FIRST-PASSAGE PERCOLATION ON LADDER-LIKE 
GRAPHS WITH HETEROGENEOUS EXPONENTIAL 

TIMES. 



HENRIK RENLUND 

Abstract. We determine the asymptotic speed of the first-passage per- 
colation process on some ladder-hke graphs (or width-2 stretches) when 
the times associated with different edges are independent and expo- 
nentially distributed but not necessarily all with the same mean. The 
method uses a particular Markov chain associated with the first-passage 
percolation process and properties of its stationary distribution. 



1. Introduction 

Consider a graph G with vertex set V and (undirected) edge set E C 
V xV. If there is an edge e = {v,v') joining vertices v and v' , they are 
said to be adjacent. A path 7r{v, v') between v and v' is an ahernating 
sequence of vertexes and edges (fo,ei,ui, . . . ,en,Vn) such that Cj = {vi^i,Vi) 
for i = 1, . . . , n, fo = u and u„ = v' . 

We think of each edge e as being associated with a (typicahy non-negative) 
random time ^e, formally we define ^ = {^e, e G E} to be a stochastic process 
indexed by the edges of the graph. We define the time T7r(v, v') of a path 
7r{v, v') to be 

Tti{v,v')= V and T{v,v') = inf Tti{v,v') 

to be the shortest time of any path between v and v' . T(v,v') is called the 
first passage time from v to v' and is the subject of investigation in first 
passage percolation, see e.g. |S W78j . 

We may also consider first passage times between sets of vertices. If Vi 
and V2 are such sets, then 

TiVuV2)= inf T{v,v') 
veVi,v'eV2 

is the first passage time from Vi to V2. 

Usually, all first passage times under consideration are from some fixed 
and given set Vq to a varying set V', in which case we denote by T{V') the 
first passage time from Vq to V' , or T{v') if V' is the singleton set {v'}, and 
refer to this simple as the first passage time of V' . 
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We can think of this as a model for a contagious disease. At time zero a 
subset Vo of the vertices are infected and subsequently the disease spreads 
to adjacent vertices. The time for the infection to be submitted from v to 
v' along the edge {v,v') is ^(^^^'). 

A typical question, and indeed the one we shall be concerned with, is how 
fast does the infection spread? This of course requires some metric on the 
set of vertices. We shall look at cases of ladder-like graphs with "height" 
as the measure of distance and independent edge times, each having an 
exponential distribution. 

Acknowledgement. The author wishes to thank Robert Parviainen for 
encouragement and help that borders on collaboration, as well as professor 
Svante Janson for interesting discussions on the mathematical technicalities 
of this manuscript. 

2. First passage percolation on ladder-like graphs 

A ladder-like graph, or width-2 stretch, is a graph whose vertex set is 
y = N X {0, 1}, N denoting {0, 1,2,.. .}, with the ordinary Euclidean metric, 
considered when each vertex is placed in the (real) plane M^. Our interest 
is on a particular class ^ of such graphs where vertices may be adjacent if 
they are no more than Euclidean distance \/2 away from each other and, in 
addition, the graph is translational invariant, see Figure [TJ By translational 
invariant we mean that we can superimpose the graph on itself by shifting 
it one unit length to the right in Figure [H 




12 3 

Figure 1. Part of the vertices of a ladder-like graph. The 
possible edges of a graph in ^ are dashed. Note that it must 
also be translational invariant. 



A vertex (x, y) in such a graph is said to be at height x S N and level 
y G {0,1}. We will consider Vq = {(0, 0), (0, 1)} to be infected at time 
zero, and we consider first passage times of the height n, i.e. the set Vn = 
{(n,0),(n,l)}. 

It should be noted that the Euclidean metric only serves as to define the 
class of graphs we are interested in. Once this is done, we only consider 
"height" as a measure of length in the graph. In this sense, the "length" 
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of the (possible) diagonal edges are not longer that that of the (possible) 
horizontal edges. 

Typically the focus is on the a.s. limit 

(7 = (7(G,4) = lim , 

n— >oo n 

called the time constant. Now, the first passage times on these graphs 
are subadditive and the almost sure convergence of T{Vn)/n relies on the 
subadditiv ergodic theorem. See Proposition 4.1 of |Sch09j for the details. 
In our calculations it is rather the asymptotic speed ^ of percolation that 
appears naturally, but this is just the inverse of the time constant, = l/cr. 

The investigation of the rate, or speed, of percolation on ladder-like graphs 
began with |FGS06] . which consider the ladder, i.e. where all vertices at dis- 
tance 1 are connected by an edge. That paper gives a method of calculating 
a when ^e^s are independent with the same discrete distribution, as well as 
a method for getting arbitrarily good bounds for the same quantity when 
the distribution is continuous (and well behaved). 

A few years later |Sch09j and [RenlOj independently studied the case of 
having independent and identical exponential times. In |Sch09j . six graphs in 

are considered and explicit expressions for a are found for three of them, 
using recursive distributional equations. It should be noted that these six 
graphs are in fact "all" graphs in the class ^ as far as a is concerned, in 
the sense that any other graph is either trivial or it has the same a as one 
of these six. 

Another method for finding a is employed in [RenlOj . which studies the 
ladder. There a Markov chain is found and the speed of percolation can be 
calculated "at stationarity" , a method that will be repeated below in a more 
general context. The method also allows the study of another aspect of the 
percolation process called the residual times, defined in |RenlO] . 

Now, we fix some notation that will be used throughout this paper. Let 

Nt = sup {x : T[{x,i)] < t} (1) 

«6{0,1} 

denote - in the interpretation of the model of a contagious disease - the 
height of the infection at time t. 
Let 

Mt = sup{x : max{r[(x, 0)],T[{x, 1)]} < t} (2) 

denote the largest height at which both levels have an infected vertex. 
Now define the front process Ft as 

Ft = Nt- Mt. (3) 

In some situations Ft is a Markov chain on N. As we shall see, the problem 
of determining the percolation speed 'f is more or less equivalent to finding 
the stationary distribution of this Markov chain. 
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We will denote by Qn ~ "i^n the property that 



lim — — = 1 

n— !>oo 



for real sequences Qn and Let z- and denote falling and rising facto- 
rials, respectively, which may be expressed as 

rfz + 1) 



z^ = z z-l) 



= z(z + 1) 



z-n + 1 - — — : 

L [z + 1 — n 

( , IN r(z + n) 
T{z) 



■, and 



3. The ladder with heterogeneous exponential times 



In this section we generalize the result of |RenlOj to a ladder with het- 
erogeneous exponential times. Let each horizontal edge be associated with 
an exponential time with intensity Ah > (i.e. an exponential with expec- 
tation 1/Ah) and each vertical edge be associated with an exponential time 
with intensity Ay > 0. Random variables associated to different edges are 
independent. Notice, since we may equally well measure intensity in units 
of Ah, there is no loss of generality in saying that Ah = 1. Let us do so and 
denote Av simply by A. 

Note that |RenlO| is concerned with the case A = 1. 

1 * * m m (ft <■» A n--- 
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Figure 2. Infected nodes at time t marked as black. Here 
Mf = 2>,Nt = & and Ft = 2>. Not shown in the picture is that 
each vertical edge has intensity A and each horizontal edge 
has intensity 1. 



Now, Ft^ defined by ([T]), ([2]) and ([3]), is a continuous time Markov chain 
on N. This fact becomes clear as we write down the intensity matrix Q of 
the process. Denote by t' = inf{t > to : Ft ^ Ft^} the first time, after to, 
that the F-process changes its value. 

Consider first the case of Ft^ = 0, which means that Mt^ and A^t^ both 
equal, say, n. Then both (n, 0) and (n, 1) are infected, and there is no 
infected node at any greater height. The only states that can be infected at 
t' are {n+1, 0) or {n + 1, 1), both resulting in Ft' = 1, and both as a result of 
the infection spreading along a horizontal edge. Informally, we will say that 
these edges lead to the state Ft/ = 1, which is rather inaccurate since edges 
connect different vertices and not different states of the front process, but a 
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convenient terminology. So, thus far we have estabhshed that the intensity 
from state to state 1 is 2 and the intensity away from state is -2. 

If FtQ =3, as in Figure [21 there are 3 vertical edges leading to states 2, 1 
and respectively. There are 2 horizontal edges leading to states 2 and 4, 
respectively. Thus, the intensities to and 1 are both A, to 2 it is 1 + A and 
to 4 the intensity is 1. 

In a similar manner, one derives the following form of the intensity matrix: 



/ -2 2 \ 

1 + A-2-A 1 

A 1-hA -2-2A 1 ... 

A A 1 + A -2-3A 1 

A A A 1-fA -2-4A 

A A A A 1+A 

V ■■• / 



This Markov chain is irreducible, so if we can find a stationary distribution 
n = (ttq, VTi, . . .), on N, such that HQ = 0, then this is the unique stationary 
distribution of Ft, and as t — )• oo, Ft will converge to it. 

Now, our main question is: what is the asymptotic speed of percolation? 
When Ftg = the intensity by which Nt increases is 2, since there are two 
horizontal edges from the state Ft^ = to the state Nf = A't^ + 1 (and no 
other edge makes any difference to either of these processes). When Ft^ > 0, 
there is one horizontal edge leading to Nf = Nt^ + 1 (and no other edge 
that makes any difference to the A^t-process) . Hence, at stationarity we get 
the speed of percolation as 

r = 2^o + l(l-vro) = l + 7ro. (5) 

Next, we proceed to find 11 through the equation HQ = 0. One way of 
doing this is to express each -Kn, n > 0, in terms of ttq; 7r„ = a„7ro — b^- 
Now, the first entry of = HQ is 

= —2^0 + TTi + A TTj. 

As J2j>o '^3 ~ ^-i which implies Ylij>i "^j = 1 ~ ^O) we get 

= - (2 + A)7ro + A + TTi yielding vri = (2 + A)7ro - A. (6) 

Repeating this procedure for the equations corresponding to columns 2 and 
3 of Q, yields, after isolating 1^2 and vra, 

= (2A2 + 7A + 2)7ro- (2A2 + 3A), and (7) 

vrg = (GA^ + 26A2 + 22A + 2)ttq - (GA^ + UA^ + 6A). (8) 

Next, taking the difference of equations corresponding to columns k and 
k + 1 for any k > 2 yields 

= TTk-i - [Xk + 3)7rfc + {\{k + 2) + 3)7rfc+i - ■Kk+2: 
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which we shah rewrite as 

-Kn = (An + 3)7r„_i - (A(n - 2) + 3)7r„_2 + vr„_3, n > 4. (9) 

We formahy define a„ and hn by the relation 7r„ = a„7ro— 6^ for n > 1. Thus, 
for n = 1,2,3 they are defined by equations ([6]), ([7]) and (l8|), respectively, 
and for n > 4 can be found iteratively through ([9]), with or hn in place 

of -Kn- 

As necessarily vrn ^ as n — ?• oo, we can find ttq as 

ttq = hm — . 

n— !>oo Qn 

This far we have copied the procedure of [RenlOj . which proceeds to prove 
the link to Bessel functions, along the lines of how this fact was originally 
discovered. Here, we start instead with the Bessel functions. 

3.1. The Bessel functions. The Bessel functions of the first and second 
kind, Ju{z) and Yi,{z), as described in e.g. jWat22| . both satisfy the recursion 

C,+i{z) + Ci{z) = — -C,{z). (10) 

z 

Define, for any (real numbers) n, A and B ^ 0, 

Jn = Jn{A, B) = Jn+A/Bi'i/B) and Yn = Yn{A, B) = Yn+A/B{'i/B). 
Then J„ and Y^ both satisfy 

C„+i + Cn-i = {A + Bn)Cn. (11) 

Define the function 

Tn = T{n,m,A,B)=^[JnYrr,- JmYn]- (12) 

This function inherits recursion (|11|) (in parameter n) 

T„+i + T„_i = {A + Bn)Tn. (13) 

Next, define 

A„ = A(n, m) = A(n, m. A, B) = T„ - T„_i. 

When m and/or A and B are clear from the context, or unimportant, 
we write A„ or /\{n,m) instead of the lengthy A(n, m, A, i?). As the next 
claim makes clear, this function is designed with ([9]) in mind. 

Claim 1. A„ satisfies 

^n = [l + A + B{n- l)]A„„i -[l + A + B{n- 3)]A„_2 + A„_3- (14) 
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Proof. By using the definition = — T^-i on the right hand side of 
(jl4p . rearranging terms so that relation (I13p can be apphed (as indicated in 
(fTSl) ). we see that the right hand side of (fT4l) equals 

[A + B{n - l)]T„,_i -2[A + B{n - 2)]T„_2 +[^ + B{n - 3)]T„_3 + (15) 
. , . ' 

+ Tn-1 — 2T.„_2 + 2T„_3 — T„_4 

= T„ - T„_i, (16) 

which is the same as the left-hand side of ()14p . □ 
Claim 2. 

T(n,n,yl,S) = and (17) 

T(n-h l,n,A,5) = 5. (18) 

Proof. The first assertion follows trivially from the definition (jl2p . 

For the second assertion, we need two important relations for the Bessel 
functions. The first is that 

zJl{z) = vJy{z) - zJu+i{z), Vv. (19) 

and the same holds with Yi, replacing J^, see 3.2(4), p. 45 and 3.56(4), p. 
66 of |Wat22j . resp. The second is a fact relating to the Wronskian, namely 
that 

— = J,{z)Yliz)-Jliz)Y,iz), \fu,z^O, (20) 

TTZ 

see 3.63(1), p. 76 of |Wat22j . 
(fT9|) inserted in ([20]) gives 



^ = tt[J,+i{z)Y,{z) - Mz)Y,+i{z)] 

and z = 2/B gives the desired result. □ 

Remark 1. It should be noted that the T -function does not solve a recursion 
of the form ( [j^[ ) with arbitrary initial conditions. Indeed, by relation il3\) 
we get 

An+i - An = iA + Bn- 2)T(n, m. A, B). 

Hence, with A and B considered fixed, two points of the sequence of A.-s 
(such that A + Bn — 2 ^ is enough to determine m - which may not be 
unique, but in general not arbitrary - which in turn determines the entire 
sequence. 



From the preceding claim, and relation ()13p . it is straightforward to cal- 
culate: 



Claim 3. 



A(m,m,^,B) = B, 
A(m + l,m, A, B) = B, and 
A(m + 2, m, A, B) = [A + B{m + 1) - l]B. 



8 HENRIK RENLUND 

Next, we need some asymptotic properties of T„. 
Claim 4. As n — )• oo, 

T(n, m, A, B) ~ J^r(n + A/ B)B'''^^'^ . (21) 

Proof. |AS64j lists the following asymptotic relations, for fixed x and n 
tending to infinity, 

V27rn V2n/ \ irn \ex J 

hence, by the well-known Stirling formula r(n + 1) ~ \/27rn""'"^/^e~" we 
conclude 

~ r{n + A/B + l)B^+^/B ^^^^ 

^ _lr(n + 
vr 

from which the claim follows. □ 

Now, we have the tools to continue with the problem of finding the speed 
of the percolation process. 

3.2. Finding the stationary distribution. Now, (fTi|) fits ([9]) if we set 

A = 2 + X and B = X. One way of trying to describe a„ and 6„ is through 
some linear combination of A:s. By inspection of Claim [3l we choose to 
work with A(n, 1) and A(n, 2). Then one finds that 

2A2 + 8A + 5 ^ , , A + 3 ^ , , 

On = 7 • A(n, 1) — • A[n, 2), and 

A A 

2A2 + 4A + 1 , , , A + 1 , , , 
bn = ^ A(n,l)-^.A(n,2), 

is true for n = 1, 2, 3 which makes it true for every n. 
Next, with the aid of Claim U we calculate 

^ = li^ ^ = (2A^ + 4A + 1) J2+2/a(2/A) - (A + 1) J3+2/a(2/A) 

n^a„ (2A2+8A + 5)J2+2/a(2/A)-(A + 3)J3+2/a(2/A)' ^ ^ 

and thereafter, a long, but straightforward, calculatiorQ yields, for n > 1, 

2 

TTn = c(Jn-l — Jn)-, with C = 

^ ' (2A2+8A + 5)Ji-(A + 3)J2 

after which it can be verified that IT is indeed a distribution since, 

n 

y^TTfc = 1 - cJn 1, 
fc=0 

as J„ ^ by (|22D . 



^The analogous calculation to that of (2.19) in [Ren 10] . 
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The speed Y of percolation in this model as given by ([5]) is 1 + ttq, with 
ttq as in (j23|) , also depicted in Figure [3] with A ranging from close to zerc0 
to 20. 

2 



1,8 



1,6 



1,4 



1,2- 



1 

5 10 15 20 

Figure 3. The speed of percolation on the ladder as a func- 
tion of (exponential) intensity A on the vertical edges (as 
measured in units of the intensity of the horizontal edges), 
plotted by Maple. Note that the speed at A = is 1. 

Although not defined, the speed as a function of A = is 1, since A = 
must be thought of as not having any vertical edges at all (or equivalent that 
the times associated with these are infinite) and as such the time to reach 
Vn is the time it takes until the first of two independent Poisson processes 
reaches n. Although this is always a bit faster than the time for a single 
such process, the effect wears off as re tends to infinity. 

Similarly, the speed at A = oo is 2, if we interpret A = oo as having 
zero time associated with vertical edges. Then infection between adjacent 
vertices on different levels is immediate, and the time it takes to move one 
step up the ladder is the minimum of two exponential variables each having 
one unit of intensity, i.e. an exponential random time having intensity 2. 

Having the intensity 2 associated with an edge is equivalent to having two 
(independent) edges associated with intensity 1, since the minimum time of 
two intensity 1 exponential random variables is an intensity 2 exponential 
random variable. Some examples of the speed of percolation on different 

*^It appears to be difficult to evaluate the speed for very small values of A, as B = 1/A 
becomes very large. In Figure O it looks as if the speed at A = is ~ 1.2. However, this 
is not the case. 
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/ N l I 7J5(2)-2J4(2) 
ya.) ± -r- 15J5(2)-4J4(2) 



(U\ 1 I 17J4(1)-3J3(1) („\ 9,0 7J7(4)-3J6(4) 

'^'-'^ 29J4{1)-5J3(2) {<^)'^-l-^ 15J7(4)-4J6(4) 



Figure 4. The speed of percolation when all edges are asso- 
ciated with random variables having intensity 1. Numerical 
values to two digits are, from left to right, 1.47, 1.59 and 
2.74. 



ladders are shown in Figure HI Notice that in the rightmost graph, we 
calculate the speed as two times that of having A = 1/2. 



4. Adding diagonals to the ladder 

Next, we see what happens when diagonals are added to the ladder. For 
this model we cannot find a solution through Bessel function, as was the case 
with the ladder. Neither can we consider arbitrary intensities associated to 
the diagonals, these must be of the same intensity as the horizontal, else the 
front process becomes intractable. 



V 




v' 



Figure 5 . Intensities associated with vertical edges have in- 
tensity A, as measured in the unit of intensity that is associ- 
ated with both the horizontal and diagonal edges. Note that 
the graph is the same if we interchange v and v' (without 
breaking any edges). 



Part of the graph we are now considering is depicted in Figure O In this 
figure we see two nodes v and v' . If we interchange the positions of these 
two nodes, without breaking any edges, the graph is effectively unaltered. 
The vertical edges are still the same (they are now "upside down", but this 
is irrelevant as edges are undirected). Four horizontal edges have become 
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diagonal and vice versa, but as these are all associated with the same in- 
tensities, nothing has essentially changed. This is, of course, the reason we 
need to have the same intensities on horizontal and diagonal edges. 

Consider the front process, defined by ([l]),([2]) and ([3]), at time t. If we 
disregard all nodes below height Mt, there are 2^ different sets of infected 
nodes (at t) that could yield {Ft = k}. From any such state the intensity to 
move to to any other state that yields {Ff = k'} is the same, and as such 
the front process will be a Markov chain on N. 

In Figure El we illustrate this when the front process is found to be 2 at t. 




Figure 6. Two of the four possible sets of infected nodes at 
t (disregarding nodes below height Mt) that would yield Ft = 
2. The infected nodes are marked as black. The remaining 
two sets are found by mirroring the cases above through a 
horizontal line. 



Now, it is straightforward to compute the intensity matrix Q of the front 
process, e.g. when in state 2, as in Figure [6l then intensities to state is 
1 + A, to state 1 is 3 + A and to state 3 is 2, and by similar considerations 
we get 



Q 



( -4 
2 + A 
1 + A 
1 + A 
1 + A 
1 + A 

V 



4 

-4- A 
3 + A 
2 + A 
2 + A 
2 + A 




2 

-6-2A 
3 + A 
2 + A 
2 + A 





2 

?-3A 
3 + A 
2 + A 






2 

-10-4A 
3 + A 



We proceed to find the stationary distribution 11 = (7ro,7ri, 
gously to the computation of ([6]), d?]), dH) and ([9]) we get 



(24) 



7r„ 



Analo- 



n > 4, (25) 



[(2 + A)n + 3]7r„_i - [(2 + X){n - 2) + 4]7r„_2 + 27r„_3, 

and if we let 7r„ = c„7ro — dn, for n > 1, we have 

(ci,c2,C3) = (5 + A,28 + 17A + 2A^226 + 226A + 68A2 + 6A^), and 
{di,d2, ds) = (1 + A, 8 + 9A + 2A^ 66 + 98A + UX^ + 6A^), 
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and Cn,dn for n > A can be determined via, relation (|25p with Cn or dn in 
place of TTn- As in the previous section, we aim to determine ttq via 

hm ^. (26) 

n— 5>oo Cn 

Once, this is done, the speed of percolation Y is 

r = 47ro + 2(l-7ro) = 2(l + 7ro). 

4.1. Turning to generating functions. Next, we employ the recursive 
techniques of [JanlO] to get an explicit expression for Cn and dn- For ease 
of comparison, we will adopt the notation of that paper. [JanlOj considers 
recursively defined sequences of the form 

K 

an = '^[ai{n - i) + Pi\an-i, n>K, (27) 
1=1 

where ai, . . . , ax and /3i, . . . , Pk are given numbers, and the initial part of 
the sequence oq, • • • , clk-i is fixed. Also, the numbers have been normalized 
as to give ai = 1. We may work in this setting if we, for n > set 

a„ = c„+i(2 + A)"" and = d„+i(2 + A)"". 

Our main interest is lim^ dnj Cn which will equal limy^ ct^jcLn- 

So, we aim to find a formula for a„ given by (p7|) where K = 3 and ai = 1, 
as = and a2 typically is negative. Set a = —a^- Then the parameters 

1 

a 



/3i = 2 + 



2 + A' 

3 7 + 2A 



2+A 2+A 



- "2TA " (2TaF " " (^TXF '"'"'^ 



(2 + A)3' 

are the same for the sequences a„ and a* , which then "only" differ in initial 
values 

, 28 + 17A + 2A2 226 + 226A + 68A2 + 6A\ 

(01,02,03) = (5 + A, — , ), and 

2 + v (2 + A)^ 

, , , ,^ . 8 + 9A + 2A^ 66 + 98A + 44A^ + 6A^ 
(ai, 02, 03) = (1 + A, 2 + ^ ' J2TW ^' 

Following [JanlO| we define the generating function 

00 

A{z) = akz^ 

k=0 
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for the sequence a„. Next, relation (l27jl . with K = 3 and parameters as 
above, yields 

A{z)q{z) = A'{z)z^p{z) + r{z), (28) 

with 

3 

p{z) = ^^QjZ*""^ = 1 — az, 

i=l 

3 



l{z) = l-^piz' = l-hz- hz'^ - /93^^ and 



i=l 

2 

r{z) = ^R^z\ with (29) 

i?o = ao, i?i = ai - /3iao, and i?2 = 02 - (1 + /3i)ai - /32a2. 

Let / C M be such that G I and p(z) 7^ on /. Theorem 3.2 (i) and (iv) 
of |JanlOj reveals that there exists a solution A to (j28p in / and any such 
solution satisfies the asymptotic expansion A{z) = J2n=o^nZ^ + 0{z^~^^), 
as z — >• 0. 

Now, of the three polynomials p, q and r it is only the latter that differs be- 
tween an and a* , so let us denote these specific polynomials r{z) = ^ RiZ^~^ 
and r*{z) = respectively. The coefficients simplify to 

i?o = 5 + A, 

^ 28 + 17A + _ 7 + Z> ,) , 7 ^ 

2+A 2+A^ ^ 2+A 

~ 226 + 226A + 68A2 -h 6A3 / 226 + 226A + 68A2 6A3 \ 

¥TW ¥TW J 

+ (2W^' + '^=(^W='"' 

iiS = 1 + A, 

_ 8 + 9A + 2A2 7 + 2A,^ 1 

R*, = -(1 + A) = = a, and 



_ 66 + 98A + 44A^ + 6A^ / 66 + 98A + 44A^ + 6A 

{2TW V ¥TW 



1,3 



6 + A 

(YTa? 

Still using the notation of jJanlOj . we implicitly define 7 and the rational 
function g{z) via 

g(^) -2 , -1 , ^ N 
-—— = z +-fz +g{z), 

z'^p{z) 
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SO that 



~ 6 + 2A , . X a{a- Pi)- h- 

7 = "-Pi = — ^—-r and g{z) = . 

2 + A 1 — az 

We fix an antiderivative G to (7 as 

G{z) = Cz + D\ii{l- az), 

where 

C = — = and D = ^ + — + /3i - a = a(2 + A) = 1. 

Now, we can write down what in [JanlOj is called the principal solution to 
(j28|) . which we consider for z < 0, 



dt 



,^(1-Qz)^(l-Zi)^ 



(1- (a + t)z)^ 



Ctz^ 
1- zt 



1- zt 



1 - {a + t)z 



i=0 



dt 



(31) 



where we by three Taylor expansions and the relation {—x)— = (— get, 
as z 0, 



2 r 

~n Jo 



i=0 

N 

E 



N 



{-a)WU 



j=0 



N 



{-a - t)"\-D - l)^z"' ^ {-Ctz^)" „i 



ml 



E 

A:=0 



fc!(l -Zt)^ 



e z 



m=0 



dt + 0(z^+^) 



i=0 j=0 fc=0 m=0 

r-oo 



m! 



poo 

/ (1 - zt)^+^+^^^-'=t'=(a + t)™e-* dt (32) 
Jo 



The integral denoted / in (j32p . equals, by Taylor expansion, as z 0, 

/■oo ^ 

/ 

1=0 







n 



N 

E 



(A; + i-7-D- 1)' ; 



Jo t'o 



m 



n 



ra"^-"e-*dt + 0(z^+^) 
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E + (^)a-«(fc + l + n)\ + 0{z^^^). (33) 



1=0 n=0 ^ ^ 

So, from ()32p and (I33p . we get the following formula 



A'' / T-v , 1 Nrri" 



^ ^^^^ E a'"""(A; + / + n)!^*+i+2fc+'+- + 0(z^+i), (34) 



m=0 n=0 

which simplifies somewhat since D = 1 in our present application. Using 
this, and letting i?o and B\ be the nonzero summands of X]^o(~'^)''-^~/-?'' 
i.e. -Bq = 1 and -^i = and by defining 

^.(., „, = f + , + „), ,35) 

n=0 

we can write, when considering Aq{^z) as an infinite sum, 

2 1 
4=0 j=0 k,l,m>0 

Further, defining Bj = if j ^ {0, 1} and Ri = ii i ^ {0, 1, 2} makes 

3 

^o(^) = E E E BjRMk,l,rn)z'^+^''+'+^. (36) 

L=0 k,l,m>0 i+j=L 

Since Aq(z) can be written on the form (j34p this allows us, by Theorem 3.2 
of [JanlOj . to identify the coefficients of ^0(2;) to those of A{z), so from ([36j) . 

3 

"^ = E E E BjR,Fi{k,l,m), (37) 

L=0 2k+l+m=N-L i+j=L 

if we follow the convention that an empty sum is zero. 
Now, to get a handle on cat, consider 



2k+l+m=M 



m 

' m 



2A;-|-i-|-m=M ' ' n=0 



(38) 
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Claim 5. 

L,:= lim ,,f:l^L ='r ,,T../~^^'^ r 

A/^oo M1M7+1+' ^ /c!r(/c + 1 + 7 + i) Jo 

(39) 

ly/iere 7 = -7 - 3 = -A/(2 + A) G (-1,0]. 

Proof. The details are given in Appendix |Al □ 
We need to find a useful expression for the integral that appears in (|39|) . 
Claim 6. Define, for n G N, 

/(n) = / (1 - x)^+"xe"^ and J(n) = / (1 - 2;)^+"e"^ dx. (40) 







Then 

r(7 + n + 1) 



I(n) 



a" 



1(0) - gJ(0) 1 " (n + 1 - m)a" 
r(7 + l) a2 2^ r(7 + l + m) 



Proof. The details are given in Appendix [Bl (Note that the result is true 
for any a 7^ but that it is essential that 7 > — 1.) □ 

4.2. Evaluating the limit. Now we return to examination of the limits Lj 
given in ([39]) . 

By (I39D, (HDD and (HH) we have, for i G {0, 1, 2}, 

-'53--7TTT^7^-rT^4 + -T55(z), (41) 



a*r(7 + l) a*+ir(7 + l) a;*+2 



where 



(-C/«) _ -C/a _ -2a 

k=0 



kl 

5, = y "^-f/"^ = _hie-^/" = -2ae-^", and 



°° fc( -C/a)^ _ C 
kl a 

k=0 
00 k+i 

=^'^ f{k,m,i), where 

k=l m=l 

/(..^.,)= '-^'X-:j;°"'*t'V'"'' - (^^) 

k!1 (7 + 1 + m) 
This yields, by rearranging ()4ip and using C ja^ = 2, 
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Recall that a G (0, 1/2] and that 7 > — 1 so that / defined in (|l2]) satisfies 
\f{k, I, m)\ < \k + i + 1 — m\/k\m\ when fc > and m > 1, and so the sum 
S{i) is absolutely convergent and we are thus free to interchange the order 
of summation 

00 fc 00 fc 00 00 

•^(0) = E E 0) = E E 0) = E E 0) 

A;=Om=l fc=l m=l m=l A;=m 

00 00 00 00 

= E E + 0) = E E + 

m=l j=0 i=0 m=l 

= y (-2ay (j + 1) y ^ (-l)'"(2a^)'" ^ ^ ^^^^ 

if we define 



^(A:) = y(-2ay(i + l + A:)y- ^ ^ -. (45) 

Similarly to how we handled ^(O) we get 

00 00 

'5(i) = EE/(^'+"^-i'"^'i) 

j=0 m=l 

'"^ i^ + l)2L.r(m + l + 7)r(m + l + j-l) 

\m (r) '2\m 



1 ^ (-ir(2a^r , 
2a {m-l)\T{m + l+^) ^' 



2a ^ (m- l)T(m + 1 + 7) 

m=l 
2(l-7)/2 

= -^^Ji+^{2V2a) + ^(1), (46) 
where J denotes the Bessel function of the first kind, as well as 

00 k+2 00 k+l 00 

^(2) = E E 2) = E E 2) + E + 2' 2) 

fc=Om=l fc=Om=l fc=0 

= EE/(^-+--i'-'2) + «^E 



A;!r(A; + 3 + 7) 

m=lj=0 fe=0 ^ " 

y 7 r + ^(2) + :^J2+y(2V2a) 

2a ^ (m-l)\T(m + l + ^) 2a^ ' 

m=l 

2(l-7)/2 ^ o-'^/^ 

2^-^ Ji+^(2V2a) + ^ J2+^(2V2a) + ^(2). (47) 
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4.3. Evaluating the limiting fraction. Now, we go back to the question 
of determining the hmit in (I26p . First, the formula for oat in ()37p equals, by 
the result in (l39l). 



UN ~ BoRoLoN''+^N\ + BiRqLq{N - lf+^{N - 1)! 

+ BoRiLiiN - iy+^{N - 1)! + BiRiLi{N - 2y+^{N - 2)1 
+ BoR2L2iN - 2y+\N - 2)! + BiR2L2{N - 3)T'+3(iV - 3)! 
~ N''+^N\BoiRoLo + RiLi + R2L2). 



Recall, from preceding equation (|35p . that Bq = \ and the Ri'.s we want to 
use, which we have denoted Ri and R^ for and a|^, respectively, are to 
be found in (j30|) . 

Note also, from (HOl), that 



1(0) + J(0) = /V - x)-^/(2+A)(^ ^ l)eX/(2+A) 
JO 



.(2 + A)(l-x)l^^/(2+A)gX/(2+A) 



x=0 



2 + A = 1/a, 



and a straightforward property of o5^(x), as defined in (l45]l . is 

y + 22 



x^(0) + 2/^(1) + zy{2) = {x + y + z)y 



X + y + z 



(48) 



Then, from dSU]), (03]), (glD, (gU]), (gT]) and (gS]), recalling also that a 
1/(2 + A), we get 



aN 



^ RiLi = (5 + A)Lo - 7aLi + Aa'^L2 



i=0 



e-^-[/(0) + J(0)] 1 



ar(7 + 1) 



+ ^ [(5 + v)S{^) - 75(1) + 4S(2)] 



-2a 



+ 



2(l-7)/2 



2-7/2 



, , . Ji+^(2\/2a) +4 ^ J2+^(2\/2a) + -^( 

r(7 + l) 2ai ' 2ai ' - ^ 



a 



a 



where we in the last equality used the result of (g8|) . Similarly we get 



''N 



iV7+iiV! 



^ R*Li = (1 + A)Lo + all + OL2 

-2a 2(^-'^)/2 1 

+ „ Ji+^(2V2a) + -^(a) 



r(7 + 1) 2aT' 



a 
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Next, we note that 

(-l)"^(2a2)™ 



fr^ ^^^r(m + i + 7)r(m + i + j) 



oo 



-2^( ^^ + ^ + "^Z^ r(m + i + 7)r(m + i+j) r(7 + i) ' 

j=0 m=U I / \ ■'J \ I I 

^ V ' 

(49) 

where we think of the sum denoted as a function of the parameter A. 




Figure 7. The speed of percolation on the ladder with di- 
agonals of (exponential) intensity A on the vertical edges (as 
measured in units of the intensity of the horizontal and di- 
agonal edges), plotted by Maple. 



Now, we have an expression for the limiting fraction 

, ^g^gy + i(^/2a)^-^J,^, + ^. 

TTn = lim = — -T , (50) 

in which we have suppressed the argument 2\f2a from the Bessel function, 
i.e. in (j50|) we should interpret as J^i^ypia). The speed 

r = 2 27rn 
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is plotted in Figure \7\ for A ranging from to 20. At zero the speed is 
3 J (V^/v^'^J {V2) ~ ^'^^ " ^^^^ exact expression is derived in Section 14.41 
below - and as A increases to infinity, the speed must reach 4. 

In Figure [8] we have given numerical values when A is 0, 1 and 2, respec- 
tively. 




Figure 8. The speed of percolation when all edges are asso- 
ciated with random variables having intensity 1. Numerical 
values to two digits are, from left to right, 2.58, 2.81 and 
2.97. 



4.4. The case A = 0. This special case has already been investigated in 
|Sch09] by other methods. In this section we show that our method gives 
the same result, since this is far from obvious. Also, the case of A = appears 
to be the only case that allows for great simplification of the expression (I50p . 
We note that for A = 0, since this implies q = 1/2 and 7 = 0, we get from 
the definition of in ([M]), 

3\ ^ (-1)^(1/2)"* 



o— n ^ ^ m—n 



i=0 ^ - m=0 ^'-i^+jy- 

fli-^y{j+l) i^yjA^) 



j=0 

„ 00 00 

= |J](-ly(^/2)V,(^/2) + J](-l)^■(^/2rl^/2iJ,(^/2) (51) 

j=0 j=0 

and since V2jJj{V2) = Jj-iiV2) + Jj+i{V2) by relation i^, we get, after 
cancellation of terms 

yo = l[MV2)-V2Ji{V2)]. 

This, together with the fact that J2(\/2) that now appears in the denomina- 
tor of ([50]) . as a consequence of relation ([TO]) , can be written as \/2Ji(\/2) — 
Jo(\/2), simplifies §Qi) to 

JoiV2)- Ji{V2)/V2 
-Jo(V2) + 3Ji(v/2)/V2' 
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SO that 

2^/2Ji(V2) 



r = 2 27ro 



3Ji(^/2)/V2- Jo(V2)' 
and thus 

1 3 



-r 4 2V2Ji(\/2) 
which is what is calculated in jSch09j . 

5. Other graphs in 

The class ^ of graphs, defined in Section [2] and Figure [H contains (at 
least) 3 more graphs of interest; these are shown in Figure [H We will 
not deal with these graphs here but make some remarks. We believe that 
the front process together with the results of |JanlOj . used in Section 
can be used to calculate the speed of percolation on the graph depicted 



in Figure 9(c) , in principle with arbitrary intensities associated with the 
vertical, diagonal and the two different horizontal edges. There would be 
added complexity as the state space of the front process would be all the 
integers, not only the non-negative ones. The speed when the intensities 
are the same is 2 tan 1-2 ~ 1-90, from the exact expression of the rate of 
percolation calculated in |Sch09j . 

The problem when trying to apply the front process to the other graphs 
is that the Markov property is lost. 




Figure 9. Three more graphs in J5f . 



Appendix A. Proof of Claim [5] 

Step 1. Let Q and K be nonnegative real numbers, a is considered fixed, 
and define 



Then we can write as 

m j j-i 
Mn = (m + i^)!^^n 



m — I 



7! m + K — i 
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We are interested in the limit of = i^^^/(m + K)\ as m tends to infinity. 
Now, as (m — i)/{m + K — i) is decreasing in i, 



n 

i=0 



m — I 



< 



m + K — i \m + K 



m 



and as a consequence, we get 



// < gOim/{m+K)) ^ 



(53) 



Next, as long as j < m, since ln[(m — i)/(m + K — i)] is a negative and 
decreasing function when < i < m, we can calculate. 



n 

i=0 



m — I 
m + K — i 



exp I ^ In 



m — I 



m + K — i 



m 



I i=0 

{m + K- 



> exp <j j In 



m — X 
m + K — X 



(m + K)"'+^ (m - j)™-* 



m 



m + K 



1 



1- i 



1-^ 



>1 



Since 



1-i 



m(m + iiT — i) ' 



the last factor in ()54p is, by Taylor expansion around 1, 



(54) 



1 - 



m{m + K — j 
and hence .s/^ is bounded below by 



+ 



K 



m{m + K — j) 



m—l 



^ ,ma/im + KW ^^ ^ ^^^^^^^ ^ + 0{l/m)). (55) 



In conclusion; from the upper and lower bounds, (|53p and (|55|) respectively. 



(m + K) ! exp 



am 
m + K 



(l + 0(l/m)). 



(56) 



FIRST-PASSAGE PERCOLATION ON LADDER-LIKE GRAPHS 23 

Step 2. Let us examine, for / > and K > —1, the fraction r(/ + K + 
l)/r(/ -|- 1) for large /. By Stirling's formula 

ni + K + 1) _ V2^{1 + KY+^+^'e-^-^ , 
r(/ + l) V2^li^y^e-^ {l + Oi[l + K] )) 

= (/ + i^y+^^e-^ ( ^ + T ) + ^^[^ + ^T'^^ 

= + + 0(1/0) (57) 

Step 3. Define the function 

l+m=M~2k ^ ' 

where £^m{k + /) is defined in (|52p . With SSi{k) defined in this way we can 
write Fi{M) from ([38]) as 

Now, examine ^i(fc) scaled by M''+'^+*+1(M - fc)!, 

■:^»(A:) _ ^ T{k + i + ^ + 1 + 1) m + l s^rn{k + l) 

jVffc+7+*+i(M- fc)! ~ r(/ -M) 71^^+0-+* ' M ' (M-/t)! ■ 

a * c 

From (j57p we get, using I + m = M — 2k, 

m^+i-n] (1+0(1/0) 

'-m) (1 + 0(1/0). 

Obviously 5 = m/M + 1/M. From (|56p we get 
am 



c = exp 



M -k 
So that, as M — )• oo. 



(1 + 0(l/m)) = exp {^} (1 + 0{l/m)). 



m=0 

-I 

[I - x)^+'+^xe"^ da;. 



And hence, since also M'=+'^+*+1(M - A;)! ~ M'^+^+^M!, we have shown the 
result in 



24 henrik renlund 

Appendix B. Proof of Claim [6] 
Recall that we have defined, for n G N, 



I{n) 



^1 -x)^+"xe"^'dx and J( 



n 



;i_2;)T'+^e"^dx, (58) 



where 7 G (—1,0] and a G (0, 1/2], although in this section it is only essential 
that 7 > -1. 

Then /(O) < 00 and J(0) < 00, since 7 > —1. For n > 1, partial 
integration reveals the recursive relationship 

-I x = l 



I{n) 



a 



x=0 



-(7 + n)(l - x) 



7+n.— 1 



X 



+ (1 - 2;)^+"]e"^ dx 







= -I{n- 1) - -J(n), 

a a 

from which follows by iteration that 



.,„).M.,o)-iE^^o-). 



(59) 



Note that the above formula holds for any n > if we interpret sums of the 
form ^™ to be zero if m < n. 

Next, we turn our attention to J(n). For ?7, > 1 we get, again by partial 
integration. 



J{n) 



n x=l 



(1 - 2;)^+'^-e"^ 
a 



J{n- 1) , 

a a 



+ 



x=0 



7 + n 
a Jo 



[(1 - x)T'+'^-ie'"^ dx 



from which follows by iteration that 



Jin) 



(7 + nY- 



a" 



^ n—l 



(7 + n)^ 



fc=0 



a" 



(60) 



Using deOj) in ((59]) yields 



.(„).(l±?li/(0) 



J(0) 



a" 



a 



,n+l 



•Si H — ^5*2 



where 



" ^-i(^ + n)!izi(^ + j)i 



nr(7 + n + 1) 



m+i) 



and 



S2 = EE 



a 



k-j 
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The above simplification of Si is easily seen to be true by the relation x- 
T{x + 1)/T{x + 1 — m). The same relation can be applied to S2 to give 

2^2^ r(7-h 1 + 7 

7=1 A:=0 ^ ' 



S2 



a" 



k) 



r(7 + n+l) ^ (n + 1 



a" 



m=l 



where the last equality follows from noting that k and j only appear in the 
summation as m = j — k, and that the values of m, i.e. 1, 2, . . ., n, appear 
n, n — 1, . . ., 1 times, respectively. 
So, we have the sought formula 

r(7 + n+ 1) 



I{n) 



+ 



a"r(7-M) 
r{j + n + l 



a 



.n+2 



r(7 + n + l) 



a" 



m 

n 

■E 

m=l 

7(0) 



n , , 

- -^0 

a 

(n + 1- 



r(7 + 1 + m) 



(n + 1 



m+i) 



^ J(0) ^^sp 

Q,2 + 1 + m) 



(61) 
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